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Guide for Running AgMIP Climate Scenario Generation Tools With R 



Version 2.3 

*,2013 



Guide for Running AgMIP Climate 
Scenario Generation Tools with R 


in Windows 


This Guide explains haw to create climate series and climate change scenarios by using the 
AgMIP Climate team ‘s methodology as outlined in the AgMIP Guide for Regional Integrated 
Assessment: Handbook of Methods and Procedures ( available for download a t mew.aKmw.one) 

Details how to: 

* install R and the required packages to run the AgMIP Climate Scenario Generation scripts 

• create climate scenarios from CMIP5 GCMs using a 30-year baseline daily weather dataset 

The Guide also outlines a workflow that can be modified for application to your own climate 
data. 


This Guide and workflow are under continuous development and we expect to add new 
capabilities and sections in future versions. The methods and scaiarios here are intended for 
AgMIP research and may not be suitable for other applications. Please contact Nicholas Hudson 
falh2106<acohunbtitedul if you have any questions or comments. 


By Nicholas Hudson fnlh2106fg)colpmbU.cdul CCSR | Columbia University 

and Alex Ruane fataanderteJnmc@iiMa.aiTl. NASA Goddard Institute for Space Studies 

With special thanks to Yunchul Jung (Univeraity of Florida) and Sonali McDennid (NASA 

GISS) 
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★Background 

Version 2.3 of the Guide for Running AgMIP Climate Scenario Generation Tools with R 
provides updated links for downloading the input data and AgMIP Climate Scenario Generation 
scripts for R. Version 2.3 also corresponds with the updated version of the script 
“agmipfarmclimate.R” that was available with version 2.1. 

The updated version of “agmip farmclimate.R” includes a correction to handle negative vapor 
pressure (Vprs) values and NaN dew point temperature (Dewp) values as well as to correct NaNs 
created when calculating the change factor for precipitation (Rain). If you have encountered any 
of these errors when running “ run agmip farmclimate.R” . we recommend that you rerun 
“run agmip farmclimate.R” with the updated version of “agmip_farmclimate.R”. 

If you have not encountered these errors you do not need to rerun “ run agmip farmclimate.R ”. 
but we highly recommend that you use the updated version of “agmip farmclimate.R” for future 
analysis. 

The updated version of “agmip farmclimate.R” can be downloaded from the links provided in 
the section Downloading R scripts and data from WebDrive . 

New sections or sections that have been substantially updated are marked with green stars (★). 
These stars have not been updated and reflect changes that were made for version 2.1. 

The only difference of version 2.3 from versions 2.1 and 2.2 is that version 2.3 contains updated 
WebDrive links for downloading the necessary scripts and minor adjustments to 
“agmip_simple_delta.R” and “agmip_simple_mandv.R” to ensure that Tmax > Tmin in the 
generated climate change scenarios. 
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Installation 

This section explains how to install R and the packages necessary to run the climate scripts. 

R is a free environment that can be used for statistical analysis and graphic development. To run 
the scripts properly, packages are required to add functionality to R. These packages are free, 
easy to download, and importantly must be loaded, but not reinstalled, every time you begin a 
new R session. The AgMIP Climate Scenario Generation scripts that are described in this Guide 
have been designed to load the required packages automatically as long as the packages have 
been properly installed. 
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Download and Install R 

Download the most recent version of R from <httn^/cran.r-profecLorg/bin/wladowi/bagefi >. 

The image below shows a screenshot from the provided link to the R website. As you can see in 
the below image, the most recent version of R is 3.0.1 (as of the publication of this Guide). If 
there is a more recent vasion of R, you should download that version. To download the 
installat ion file in your web browser, click on the link indicated by the red arrow. Run the 
downloaded .exe file and follow the steps of the installation wizard. After the installation is 
complete, open R. 



tttft teatufts f) Ete Yereop 

If you want to double-check that the package you have downloaded exactly matches the package distributed by R. you can compare the mdSsum of 
the exe to the true fmaenamt You will need a % eruon of md5*um for window 1 * both graphical and command line vers toes are available 

Frequently asked question* 

• Berio ij» «n R ufrro Mini Wnvfoim Yiiq- 

• H e" j 9 1 m my KQYloai y £ T?.«° pf .V 

• Should 1 run 32-bti or 54-bti R'* 

Please see die R FAQ for general information about R and the R Windows FAQ for Windows -specific information 

Other builds 

• Patches to this release are incorporated in the t-paichcd snapshot build 

• A build of the development muon (winch will eventually become the next major release of R) » a\ ailahle m the r -dev el snapshot build 

• Ptei ioUf 

Nose to webmasters A stable link which will redirect to the current Window*! binary release tt 
<CRAN MIRROR-"* bin windowx base release him 


Last change 2013-05- Id. by Duncan Murdoch 


Note: The R website provides a lot of helpful information if you are experiencing any 
difficulties downloading or installing R onto your computer. Also, there are versions of R 
available for Mac OS < httn://cran.r-pro1ectorg/bln/mmco»x/> and Linux < http://cran.r- 
Drolert.org/hln/HniiTfa ’ Again, the R website is very useful so if you are installing R on 
either of these platforms and are experiencing any difficulties, your questions may be 
answered by information provided by Ibis website. 
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Download and Install Packages 

1. After installing R onto your computer, you will need to install the “R.matlab”, “R. utils”, and 
“MASS” packages. Let’s start with, the “R.matlab” package. You can perform this action by 
either 

• Selecting R: Packages-> Install package(s). . . from the menu bar (shown by the red arrow 
below) or by 

• Typing the command 

install. packages ("R.matlab") 

into R Console (the orange arrow below). 


U RGui {64-bit} 


i CiD 
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After selecting the nearest CRAN mirror, R will install the selected package “R.matlab” as well 
as the dependent packages “R.oo” and “R.methodsS3” to your package library. Occasionally the 
closest mirror does not work. If this is the case, simply try a different location. A successful 
install will display the following message in R Console: 

package 'R.oo' successfully unpacked and MD5 sums checked 
package 'R.methodsS3' successfully unpacked and MD5 sums checked 
package 'R.matlab' successfully unpacked and MD5 sums checked 


IMPORTANT: The next step, Step 2, is crucial for the installed packages to run correctly. 
The packages must be loaded, but not reinstalled, every time you begin a new R session 
to perform the operations required by the provided climate scenario scripts. However, 
the provided scripts described in this Guide automatically load the necessary packages with 
the 3 line code shown below: 

library <- c ("R.matlab", "R. utils") 

lapply (library, require, character .only = T) 

rm (library) 

Thus, if you will only be running the provided climate scenario scripts, you can skip Step 
2 and proceed to Step 3. If you are designing your own scripts and would like to retain the 
functionality of the provided scripts you will either need to incorporate the three line code 
above or follow Step 2 to load the packages. 
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2. Load “R-matlab” package. You can perform Ihis action in one of two ways. 

• You can select R: Package s-> Load package. . . from the menu bar (shown by the red 
arrow below), which brings up a pop-up box titled “Select one.” Choose “Rjnatlab” as 
indicated by the orange arrow. 

• Alternatively, you can type the command 

library (R. mat lata) 

into R Console as is indicated by the yellow arrow below. 


R 

f*. t<M Jtww Mile IgntaQW 1 b«*p 

load pactaga 

Set CRAN nwroT. 

Select icpotdcnrt. 
hneui pa<keg*t»l- 
Update packages.. 

IniuM package* jf from local tip Met... 


■ o| B| 


U R Cental* 


Please aeleci 

also mat ailing i 


n fes 


/ vindeva/ cont rib / 3 . 0/R . oo_f 
43 Ebl 


trying URL 'http; 

Content typo ' apt 
opened URL 
downloaded *43 Kb 

trying URL 'http://wataon.bci .nin.gov/cxes_mirror/bin/vtndove/cantrib/3.0/R.B«t5 
Content type 'application/ tip' length 516*5 bytea <50 Kbl 
opened URL 
downloaded SO Kb 

trying URL 'http: //vatacn. nci.nih.gov/eran_nirror/bm/vindcva/contrib/ 3. 0/R.satf 
Content type 'application/np' length 14314*’ bytea |13» Kbt 
opened URL 
downloaded 135 Kb 

package 'R.oo* aucceaafully unpacked and M£>5 auns checked 
package 'R.sethoda33* ruccesafully unpacked and KD5 iju checked 
package 'R.natlab' aucceaafully unpacked and MDS iju checked 

The downloaded binary package a are in 

mr directory here. 

library •R.natlab *| 


mathodi 

n»g<v 

nlm* 


R.methodsS3 

R .00 

Rut* 

rpart 

rstudio 

spatial 

tpknn 

dab 


tdtk 

took 

translations 

ut<h 

I OK 1 
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When loading R.matlab, R will automatically load the dependencies, R.oo and R.methodsS3. R 
Console’s output will inform you if all three packages were successfully loaded. The red arrow 
below points to the lines indicating “R.oo” and “R.methodsS3” were successfully loaded, while 
the orange arrow points to the “R.matlab” line. 



3. Repeat steps 1 and 2 (pages 6-8) to install and load the packages “R.utils” and “MASS”. 


IMPORTANT: The “MASS” package will be used by “agmip_simple_mandv.R” and has 
been added to the list of required packages for Version 2.0 of the Guide. If you are 
intending to run the script “run agmip simple mandv.R”, you will have to install the 
“MASS” package. Otherwise it is not necessary to install this package. 
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Downloading AgMTP Climate Scenario Generation Scripts and Data - 10 

★Downloading AgMIP Climate Scenario Generation Scripts and Data 

The following section explains: 

1. The folder hierarchy necessary for c lima te scripts and data to function properly 

2. How to download the climate scripts and. data 

3. Where to place these files within your newly created folder structure 


Creating folders for R scripts and Input data 


AgMIP c lima te scripts require the below directory structure. We recommend that you use this 
folder hierarchy as this is the easiest approach for running the AgMIP Climate Scenario scripts. 
The hierarchy of the folders should Look as follows (where blue - folders, purple - files, green - 
description, red = important information): 



Location of scripts 


l ocation of data 


1 1 R scripts, 100 KB 


dim files 


Union 


| 360 Hies. 

-4.M OB 

62 files, 

746 KB 


WorldChni 


Contains WorldClim riles i.sublnt sublon. subAlt. subPrct. and 
subTnieaii) for cadi siibremon 
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We now provide you with the folder structure already created. The first link of the next section. 
Downloading R scripts and data from WebDrive . is the zipped folder structure that already 
contains the R scripts, test files, AgMIP Excel file template, WorldClim datasets and latitude and 
longitude datasets for GCMs in the correct folder locations. Additionally, you will have to 
download the GCM climate data and place these files into “~/R/data/CMTP5/climfiles”. 


IMPORTANT: R is case sensitive so it is important to preserve the downloaded folder 
hierarchy for the AgMIP Climate scripts to function correctly. 


^Downloading R scripts and data from WebDrive 

The next step will be to download the climate data and R scripts from WebDrive (links provided 
on the next couple pages). This is NASA Goddard’s Online File Depot and has been set up by 
Alex Ruane. To gain access to these files, use the username and password provided below. 


WebDrive Log In Information 

Username 

aruane 

Password 

AgMTPRRT s 

Expiration Date 

November 23, 2013 


The provided links have a 30 day expiration date and the files were posted on October 24, 2013, 
and thus will not be available after November 23, 2013. If you are attempting to download these 
files after this date, contact either Alex Ruane at alexander.c.niane@nasa.gov or Nicholas 
Hudson at nih21 06@columbia.edn . 

You will be downloading GCM files that use RCP 4.5 and RCP 8.5 as well as the latitude and 
longitude files needed to properly read the GCM files. 

Not all of the files that you will be downloading will be described in this Guide, but they will be 
discussed in future versions. 

Below are the links for the required input data and R scripts along with a brief description of the 
files and the folder destination where the files should be placed within your newly created folder 
structure: 
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Downloading AgMIP Climate Scenario Generation Scripts and Data - 12 

R scripts, test climate files, and AgMIP Excel File Template download location: 
https://webdrive.gsfc.nasa.gov/longauth/600/nicholas.hudson/HErgQoC 

1. The above link is the zipped folder structure containing the R scripts, test files, AgMIP 
Excel file template, WorldClim datasets and latitude and longitude datasets for GCMs in 
the correct folder locations. We advise that you download the most recent versions of 
these scripts to ensure optimal performance and most robust analysis. 


Script Name 

Brief description of functionality 

acrfindspot.R 

Returns the i and j coordinates of a given location given its 
latitude and longitude and a models lat and Ion arrays 

agmipfarmclimate.R 

Formerly “acr_agmip022.R” and “acr_agmip021.R”, 
produces monthly climate files from WorldClim data for n 
farms around the seed weather station. UPDATED for 
Version 2.1 

run agmip farmclimate.R 

Runs “agmip_farmclimate.R” to produce monthly climate 
files in AgMIP format. 

agmip_simple_delta.R 

Creates delta scenarios from CMIP5 GCMs and BCSD. 

run agmip simple delta.R 

Runs “agmip_simple_delta.R” to produce basic future 
scenario files in .AgMIP format. 

agmip_simple_mandv.R 

Formerly “acr agmiplOO.R”, adjusts a climate time series for 
AgMIP mean and variability scenarios. 

run agmin simple mandv.R 

Formerly “acr_agmipl20.R”, runs “agmip simple mandv.R” 
to produce climate time series in AgMIP format. 

agmip_sirnple2full.R 

Converts basic future scenarios (Srad, Tmax, Tmin, Rain) 
output from “agmip_simple_delta.R” and 
“agmip_simple_mandv.R” into full future scenarios with 
relative humidity-controlled vapor pressure based upon daily 
Tmax. 

run agmin simnle2fulI.R 

Runs “agmip_simple2full.R” to produce full future scenario 
files in .AgMIP format. 

acr_agmip2metric.R 

Produces a record of climate metrics in a given season 

run_metrics.R 

Runs “acr_agmip2metrics.R” 


The R scripts are all located in the directory “~/R/r”. This zip file also comes with a test baseline 
file, “USAMOXXX. AgMIP” located in “~/R/data/Climate/Historical” and an example future 
scenario created from “run agmip simple delta JR ” and “ run aemip simple2fuli.R ” that can 
be found in “~/R/data/Climate/Fullscenario”. 



Handbook of Climate Change and Agroecosystems Downloaded from www.worldscientific.com 
by NASA GODDARD SPACE FLIGHT CENTER LIBRARY on 04/14/15. For personal use only. 


402 


Guide for Running AgM IP Climate Scenario Generation Tools With R 


Downloading AgMIP Climate Scenario Generation Scripts and Data - 13 

Additionally, we have provided a template to create .AgMIP historical baseline files. This 
template is called “AgMIP Excel File Template v2.0.xls” and is located in the folder 
“~/R/data/Climate//Historical”. This Excel template file uses a Visual Basic macro to create 
properly formatted base files and has been designed to place the newly created base files into the 
same directory in which the Excel template file is located. Therefore if you will be using 
“AgMIP Excel File Template v2.0.xls” to create your baseline .AgMIP files, the Excel file will 
not have to be moved. Please note, to use this feature in Excel, you must run Excel as an 
“Administrator”, which can be done by right-clicking the Excel icon and selecting “Run as 
Administrator”. 

Lastly, the zip file contains WorldClim datasets in “~/R/data/WorldClim” and latitude and 
longitude datasets for GCMs in “~/R/data/CMIP5/latlon”. 


NOTE: We recommend that you keep the zip file containing the R scripts, the Excel 
template file and the .AgMIP test file as a back-up in the event that these files become 
corrupted. This is not necessary for the other zip files you will be downloading in the 
remaining four links . 


2. Place the remaining files listed below (all CMIP5 input files including meanpr, Iwetprl, 
meantasmax, meantasmin, stdtasmas, and stdtasmin) in “~/R/data/CMIP5/climfiles” 

Mean precipitation and the frequency of wet days (meanpr and fwetprl) file download location: 
https://webdrive.gsfc.nasa.gov/longauth/600/nicholas.hudsnn/vXZ7.iiYr 

Mean maximum and minimum temperatures (meantasmax and meantasmin) file download 
location: 

https://webdrive.gsfc.nasa.gov/longauth/600/nicholas.hudson/0erCIC3 

Standard deviation of maximum and mini mum temperatures (stdtasmax and stdtasmin) file 
download location: 

https://webdrive.gsfc.nasa.gov/lonEauth/600/nicholas.hudson/xdIPIC3 

The files contained in these zipped folders should all be placed in “~/R/data/CMIP5/climfiles” 
such that there are no folders located in “~/R/data/CMIP5/climfiles”, just the 360 climate files. 
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★Creating baseline .AgMIP files with “AgMIP Excel File Template_v2.0.xls” 

Before we begin running R, we will walk through an explanation of how to create baseline 
.AgMIP files with “AgMIP Excel File Template_v2-0.xls”. We are currently seeking to create a 
script that will help in the historical period gap filling and bias correction according to the 
AgMIP protocols described in Section 4 of the AgMTP Guide for Regional Integrated 
Assessment: Handbook of Methods and Procedures , but it is useful to be able to create 
.AgMIP files in Excel regardless. If you prefer to create the baseline .AgMIP files manually, 
Appendix 1 provides a detailed outline for creating these files in Excel. 

The baseline .AgMIP file should contain the 31-year (1980-2010) daily historical climate data 
for a particular station or set of climate data. This section will explain how to: 

1 . Open the Excel file and enable macros 

2. Enter your climate data into the spreadsheet and create a baseline .AgMIP file using the 
macro enabled functionality 

3 . Create baseline .AgMIP files without the macro 


★Opening “AgMIP Excel File Template.xls” 

The provided baseline AgMIP creation template was designed to quickly create properly 
formatted .AgMIP files to be used throughout the remainder of this Guide. If you have already 
designed your baseline files, you can skip this subsection. 


NOTE: This template Excel file automatically produces baseline .AgMIP files within the 
folder in which “AgMIP Excel File Template_v2.0.xls” is saved. Accordingly, we 
recommend that you leave this template file in “~/R/data/Climate/Historical” so that the 
created .AgMIP baseline files are placed in the correct folder location and do not have to 
be moved later. 


The spreadsheet includes a Visual Basic macro (essentially an Excel script designed to make 
complicated or redundant tasks easier) which has been designed to assist with the creation of the 
.AgMIP file. This spreadsheet’s macro is compatible with all Windows versions of Excel 
created since 2003. If you are using a Mac and running Excel 2008, you will not be able to run 
this macro as Mac’s 2008 version of Excel does not use Visual Basic. All other versions of 
Excel should be able to run this macro. If you would like to run this macro, please open Excel as 
an “A dmini strator”, which can be done by right-clicking the Excel icon and selecting “Run as 
Administrator”. If you are unable to run the macro, or do not wish to enable this functionality, 
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the template can still be used to produce properly formatted baseline .AgMIP files as described 
in the subsection, Creating baseline .AgMIP files without the enabled macro . 

To retain the ability of the macro, you will have to enable the macro when you open “AgMIP 
Excel File Template_v2.0.xls”. Each version has a different way of enabling the macro. 
Generally, when you open the Excel file, a security alert will appear prompting you to decide 
whether to enable the macro or disable this feature. Enabling the macro will retain the 
functionality that is discussed in the next step. If you disable the macro, you can skip the next 
subsection and proceed to the subsection. Creating baseline AgMIP files without the enabled 
macro, to learn how to use the template without the macro. 

If you are having difficulties enabling the macro, the Windows Office website has helpful 
information on enabling macros for Excel versions 2003. 2007 and versions 2010 and 2013 . 


Creating baseline .AgMIP files with the enabled macro 

After enabling the macro, you should see the cover sheet, named “Input Data”, of the 
spreadsheet. Essentially all you have to do to create the baseline .AgMIP file is fill in all of the 
light green fields. Specifically: 

1. Enter the information for your baseline file header into the "User Entered Values" section 
on the left in the fight green fields. This includes filling out information for the header of 
your .AgMIP file to be created. We highly recommend that you fill in these fields to 
properly document the data as the AgMIP Climate scripts have been designed to 
automatically use this information. There are six “User Entered Values” required for the 
header which should be unique for each station’s set of climate data. 

8 digit baseline .AgMIP file name - The first 4 digits are used to describe the location of 
your climate (meteorological) station. The first 2 digits are an abbreviation for the 
country, which can be found at the website <http://www.web-l.com/coimtrv-codes/ > 
(e.g. US = United States). The second 2 digits refer to the specific site location (e.g. AM 
= Ames, Iowa). The fifth through eighth digits (OXXX) should remain unchanged. For 
example, a file created for Ames, Iowa, USA would be named “USAMOXXX. AgMIP” 
and “USAMOXXX” would be entered into this field. The naming convention will be 
described in further detail in the subsection, Explanation of ontput files . 
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Creating baseline .AgMIP files with “AgMIP Excel File Template_v2.0.xls” - 16 
Latitude (LAT) - latitude of the station in degrees North 
Longitude (LONG) - longitude of the station in degrees East 
Elevation (ELEV) - elevation of the station in meters 

Thermometer height (REFHT) - the height of the station thermometer in meters above 
the entered station elevation 

Anemometer height (WNDHT) - the height of the station anemometer in meters above 
the entered station elevation 

2. Arrange your 31-year (1980-2010) daily climate dataset in the proper column order 
(where Srad = solar radiation, Tmax = temperature maximum, Tmin = temperature 
minimum, Rain = precipitation. Wind = wind speed, Dewp = dew point, Vprs = vapor 
pressure, and Rhum = relative humidity) into the table in the bottom left comer of the 
sheet “Input Data” (in the light green fields). 

You should fill missing values with -99. If you do not have data for Wind and/or Rhum, 
you can retrieve this data from the Ag-MERRA dataset. Dewp and Vprs will require 
additional conversion and can be left as missing (entered as -99) for the purposes of the R 
routines described in this Guide. 

Upon entering the data correctly into these columns, the orange boxes’ color should turn 
light green and display “Yes”. 

3. Once the data has been arranged correctly, simply press the green button titled “Create 
.AgMIP file” to produce your new baseline .AgMIP file. This new file will be created in 
the same file you have positioned “AgMIP Excel File Template_v2.0.xls”. Confirm the 
presence of your new baseline .AgMIP file in “~/R/data/Climate/Historical”. If you have 
successfully produced your baseline .AgMIP file, you can proceed to the next section, 
Henning AgMIP Climate Scenario Generation scripts . 
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Creating baseline .AgMIP files without the enabled macro 


NOTE: This subsection is only for individuals who have not been able to enable the 
macro or who have opted not enable the macro. If you have successfully created your 
.AgMIP baseline file, you can proceed to the next section. Running AgMIP Climate 
Scenario Generation scripts . 


If the macro will not work or you have chosen not to enable the macro, the template can still 
useful for you. Follow steps 1-2 in Creating baseline .AgMIP files with the enabled macro to 
correctly enter your data into the first sheet (Input Data) of “AgMIP Excel File 
Template_v2.0.xls”. 

After you have entered your header information and climate data into the light green fields, go to 
the second sheet, Baseline Data. This shows the properly formatted data and header. All you 
have to do to create the .AgMIP file is to save this sheet. To do this, follow these steps: 

1. Select “File > Save As” 

2. Select the location “~/R/data/Climate/Historical/” 

3. In “File name:” enter the name of your baseline file in quotes. That is, enter the name 
you have entered into the field 8 digit baseline .AgMIP file name on the sheet “Input 
Data” with the .AgMIP extension. For example, if you entered “USAMOXXX” into that 
field, enter “USAMOXXX.AgMIP” into the “File nam e:” field including the quotes. 

4. In “Save as type:” select “Formatted Text (Space Delimited) 

5. Select “Save” 

6. A warning will appear alerting you that this format will not support multiple sheets . 
Select “OK”. 

7. There will be another warning about compatibility issues. Select “Yes”. 

8 . Your file should now be saved, however, you will want to close the now altered Excel 
spreadsheet and select “Don’t Save” to ensure the initial “AgMIP Excel File 
Template_v2.0.xls” is not altered in this process. 

The macro actually works around this “Don’t Save issue” by copying and pasting the data from 
the second sheet (Baseline Data) onto a new workbook. This can be difficult though as older 
versions of Excel do not have the functionality of Paste Special as the newer versions (both 
numerical and columnar formatting need to be retained). We encourage users to follow the 
above workflow and only paste the data into a new workbook if they feel comfortable ensuring 
that the formatting of the data is not lost in the pasting process. 
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Confirm the presence of your new baseline .AgMIP file in “~/R/data/Climate/Historical”. If you 
have successfully produced your baseline .AgMIP file, you can proceed to the next section. 
Running ApMTP Climate Scenario Generation scripts . 
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Running AgMIP Climate Scenario Generation scripts 
At this point, you should have: 

• Downloaded and installed R and necessary packages as explained in the section 
Installation 

• Downloaded the input data and R scripts placing these files into the proper folder structure 
as described in the section Downloading AgMIP Climate Scenario Generation Scripts 
and Data 

• Created your baseline .AgMIP riles as outlined in the section Creating baseline .ApMTP 
files with “AgMIP Excel File Template v2.0.xls” 

If you have accomplished the above, you are now ready to run the AgMIP Climate scripts in R. 

The following section will explain: 

1 . The input variables that will be modified for each of the run scripts 

2. How to modify the run scripts 

3 . How to source the run scripts 

4. An explanation of the output files 

The subsection Explanation of input variables to be changes in the run scripts details the 
input variables that must be defined in the run scripts, the subsection Modifying variables in the 
run scripts generalizes how to prepare the run script and the subsection Sourcing the run 
scripts demonstrates how to run the script If you have any questions with the naming 
convention imposed on the output files, please refer to the subsection Explanation of output 
files . 

A Explanation of input variables to be changed in the run scripts 

We will be changing the run scripts to reflect the information for your specific location and 
analytical requirements. The run scripts all begin with “run” and call the scripts that are 
identically named but do not begin with “run” (e.g. “ run HPmin farmcliinate.R ” calls 
“agmip_farmclimate.R” and “run agmip simple delta.R ” calls “agmip_simple_delta.R”). 
These scripts have been designed so that you only have to edit the run scripts and should not 
have to open the other scripts. 
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Essentially, the ran scripts set the variables and parameters of the analysis. These run scripts are 
also designed to loop through multiple scenarios, thus producing many .AgMIP files with 
multiple scenario changes imposed upon the baseline data (different emission scenarios or future 
time periods, for example). 

All of the run scripts have been similarly designed and share many of the same input variables. 
To simplify the explanation for using these AgMIP Climate scripts, we have provided a list of 
the input variables in this subsection that will need to be changed prior to running, or "sourcing”, 
of the scripts. To clarify, each run script will require that some, but not all, of the input variables 
listed in this subsection be defined. 

We recommend that you read through this subsection to become familiar with the input variables 
that you will be c hang in g . It is our intention that this subsection will be used as a reference later 
if you have a question about a particular input variable. In the next subsection, Modifvint 
variables in the ran scrinta. we will describe how to change these input variables. Below ia a 
screenshot of the “Input Variables” section of the script “run agmip simple ddtaJR” . 


R RGw. (64- bit! 


: 'mUmd ' 


U Mrirtor iuS 

IIMMIMMMIMIMIMIMMMMIMMMIMIMIMMMIMMIMMIMMIMMMIMIMIMIMMMIMMMIMt . 

II Input variables 

Ml These should be the only variables you will have co change for running the script for 

III another file as long as all files are located in the correct folders 

basefile | <- •U5AM0XXX* •• . AgHXP base file name in -\\R\ \data\ tCliaateWMlstorical 

basedecs <- e(1980, 2009) If Tune scale of basefile 

Ml rua.gcss sets the GCM scenario loop. Currently set to run all GCKs |li20) where 


III 1 • ACCE331-0, 

**» - 

ttt 9 - SGFBL-E3M2K, 

Ml 13 - IP3L-CX5A- LR, 

Ml 17 • MFX-E5M-LR, 

run. gats <- 1:20 


2 • bec-essd-1, 

€ - CtSHl-BGC. 

10 - Mad3£M2-CC, 

14 - IP3L-CKSA-MR. 
18 • MPI-ESK-KR, 


3 - BWC-EBl, 

1 • CSIRO-MW-6-0, 
11 “ HadS£M2-ES, 

15 - HI ROC 5. 

19 • KRI'CGOO. 


4 - CanZ5M2, 

9 • 0FBL-E5M2. 
12 ■ in»cB4 r 
14 - MIROC-ESK, 
20 - RorESHl-K 


Ml run. reps sets the RCP scenario loop. Currently set to run RCP 4.5 (9) and RCP 8.5 <5) where 
III 1 - historical, 2 - RCP 2.6, 3 - RCP 4.5, 4 - RCP 6.0, 5 • RCP 8.5 
run. reps <- c(3,5) 

III run. decs sets the Time scenario loop. Currently set to run all time periods (1:3) 

M* where l - Mear-ter* C2010-2039I . 2 - Kid- Century <2040-2069), 3 - End- .310-2099) 

run. decs <- 1:3 

III You must enter the location of the R folder into roocOir below using \\ between folders. 

Ml For exaxple, *Ci WUseraWYeur Naxe HereWDesbtopWRW 

roocDir <- •••• your directory here •••\\R\\* •• <- Enter location nere <- 

MMMMMIMM You should not have to adjust any of the variables below this line IIIIIMIMIMI • 
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Colored arrows identify the input variables that will need to be adjusted for this run script. The 
“#” symbol is used to “co mm ent out” lines and provides you with information about the input 
variables. The script skips lines that have been commented when it is being sourced. 

The next few pages describe the input variables that you will be required to adjust prior to 
sourcing the run scripts. 

rootDir - This variable, as identified by the blue arrow in the screenshot on the previous page, 
sets the location, or root directory, of your “~/R/” folder discussed in the subsection. Creating 
folders for R scripts and input data . If you created the recommended folder structure, defining 
this variable is straightforward. There is a line in every run script that reads: 

rootDir <- ’*** your directory here ***\\R\\’ 

Simply replace the text “*** your directory here 1 "**” with the location of your newly created R 
directory in the script using double forward slashes (\\) to separate folders. For example, if the 
location of the R folder is 

C:\Users\Your Name Here\Desktop\R\ 

then the rootDir should look as follows 

rootDir <- 1 C : \\Users\\Your Name Here\\Desktop\\R\\ ' 


IMPORTANT: The variable rootDir should not change from run to run. We recommend 
that you save the run script with updated rootDir to reflect the location of “~/R/” so you do 
not have to enter this information for every run. The other input variables should be 
changed with each new run to reflect the specific parameters of your current analysis. 


basefile - This is the 8 digit baseline .AgMIP file name described in the subsection Creating 
baseline .AgMIP files with the enabled macro and identified by the red arrow on the previous 
page. The first four digits describe the location of the station/site while the last four digits denote 
the type of scenario that was imposed on the historic data. Specifically, the first 2 digits are an 
abbreviation for the country, which can be found at the website < http ://www. web- 
l.com/countrv-codes/ >. while the second 2 digits refer to the specific site location and are up to 
you to assign. For a baseline file, the last four digits should be “OXXX”. For example, we have 
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defined the basefile for Ames, Iowa, USA to be “USAMOXXX”. Thus the run script input 
variable line would look like 

basefile <- 'USAMOXXX' 

More information about the n amin g convention can be found in the subsection. Explanation of 
output files . 

seedflle - Similar to basefile, seedfile is the 8 digit code of your baseline climate file. The input 
variable seedfile is used in the script “run agmip farmclimate.R” as this will be the dataset 
from which other local farm sites’ climate files will be constructed. For example, the 8 digit 
code of the seedfile would be “KEMBOXXX” if the baseline AgMIP file was named 
“KHMBOXXX.AgMIP”. 

shortregion - Another input variable used exclusively with “ run aymin farmclimate.R ”. 
shortregion is a 2 digit identifier used for a particular region for which you will be producing 
multiple files. Typically, it is sufficient to simply strip off the country code. For example, the 
country code “KE” can be stripped off the seedfile “KEMBOXXX” leaving the regional code 
“MB”. The resulting output files from “run aemip farmclimate.R” will be named 
“MBO 1 OXFX. AgMIP”, “MB020XFX.AgMIP”, “MB030XFX.AgMIP”... 

end.code - This input variable, used exclusively in “run aomin simnle2fiill.R ” denotes the 
last two (seventh and eighth) digits of the files to be adjusted. For example, if you want to adjust 
all of the files of the structure “USAM**XA.AgMIP” with a baseline file of 
“USAMOXXX.AgMIP” to the full format, the 

would be “XA”. This would set the run script to cycle through the scenarios (defined by , 
run. reps, and run.decs) of all the files in “~/R/data/Climate/Simplescenario” that begin with 
“USAM” and end with “XA”. 


basedecs - This variable, also identified by the red arrow in the screenshot, assigns the 
beginning and end year of the dataset’s 30-year time interval. Typically, the test file runs from 
January 1, 1980 to December 31, 2009, thus the years entered for the test file’s 30-year time 
interval are 1980 and 2009. To assign 

seedfile - Similar to basefile, seedfile is the 8 digit code of your baseline climate file. The input 
variable seedfile is used in the script “run_agmip_farmclimate.R” as this will be the dataset 
from which other local farm sites’ climate files will be constructed. For example, the 8 digit 
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code of the seedfile would be “KEMBOXXX” if the baseline AgMIP file was named 
“KEMBOXXX. AgMIP”. 

shortregion - Another input variable used exclusively with “runagmipfarmclimate.R”, 
shortregion is a 2 digit identifier used for a particular region for which you will be producing 
multiple files. Typically, it is sufficient to simply strip off the country code. For example, the 
country code “KE” can be stripped off the seedfile “KEMBOXXX” leaving the regional code 
“MB”. The resulting output files from “run agmip farmclimate.R” will be named 
‘MBO 1 OXFX. AgMIP”, “MB020XFX.AgMIP”, “MB030XFX.AgMlP”... 

end.code - This input variable, used exclusively in “run_agmip_simple2full.R”, denotes the 
last two (seventh and eighth) digits of the files to be adjusted. For example, if you want to adjust 
all of the files of the structure “USAM**XA.AgM3P” with a baseline file of 
“USAMOXXX. AgMIP” to the full format, the 

would be “XA”. This would set the run script to cycle through the scenarios (defined by , 
run.rcps, and rnn.decs) of all the files in “~/R/data/Climate/Simplescenario” that begin with 
“USAM” and end with “XA”. 


basedecs to run from 1980 to 2009, this input variable would be assigned as 

basedecs <- c (1980,2009) 

If you have data for 2010, you have the option of either selecting the time period from 1980 to 
2009 orl981to2010. If you chose to select 1 98 1 to 20 1 0, you would enter the input variable as 

basedecs <- c (1981,2010) 

to assign 

seedfile - Similar to basefile, seedfile is the 8 digit code of your baseline climate file. The input 
variable seedfile is used in the script “run agmip farmclimate.R” as this will be the dataset 
from which other local farm sites’ climate files will be constructed. For example, the 8 digit 
code of the seedfile would be “KEMBOXXX” if the baseline AgMIP file was named 
“KEMBOXXX. AgMIP”. 

shortregion - Another input variable used exclusively with “run agmip farmclimate. R”, 
shortregion is a 2 digit identifier used for a particular region for which you will be producing 
multiple files. Typically, it is sufficient to simply strip off the country code. For example, the 
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country code “KE” can be stripped off the seedfile “KEMBOXXX” leaving the regional code 
“MB”. The resulting output files from “run_agmip_farmclimate.R” will be named 
“MBO 1 OXFX. AgMIP”, “MB020XFX.AgMIP”, “MB030XFX. AgMIP”... 

end.code - This input variable, used exclusively in “run_agmip_simple2full.R”, denotes the 
last two (seventh and eighth) digits of the files to be adjusted. For example, if you want to adjust 
all of the files of the structure “USAM**XA. AgMIP” with a baseline file of 
“USAMOXXX.AgMIP” to the full format, the 

would be ‘XA’\ This would set the run script to cycle through the scenarios (defined by , 
run.rcps, and run. decs) of all the files in “~/R/data/Climate/Simplescenario” that begin with 
“USAM” and end with “XA”. 

basedecs to run over this time scale. Generally you should not have to change this input 
variable. 


headerplus - This variable allows you to enter additional information into the first line of the 
produced AgMIP files. Typically it is sufficient to simply enter in the location of the site (Ames, 
Iowa, USA or Embu, Kenya, for example). Occasionally it is useful to include information 
about where the file data is from or how the file data was calculated. You should not have to 
change this input variable, but feel free to include additional infoimation as you see appropriate. 


run.gcms - This variable, identified by the orange arrow, is used to control which GCMs will be 
used to create the output files. Below are two tables of the GCMs with their corresponding 
value. 


GCM 

run.gcms value 

ACCESS1-0 

1 

bcc-csml-1 

2 

BNU-ESM 

3 

CanESM2 

4 

CCSM4 

5 

CESM1-BGC 

6 

CSIRO-Mk3-6-0 

7 

GFDL-ESM2 

8 

GGFDL-E SM2M 

9 

HadGEM2-CC 

10 


GCM 

run.gcms value 

HadGEM2-ES 

11 

inmcm4 

12 

IPSL-CM5A-LR 

13 

IPSL-CM5A-MR 

14 

MIROC5 

15 

MIROC-ESM 

16 

MPI-ESM-LR 

17 

MPI-ESM-MR 

18 

MRI-CGCM3 

19 

NorESMl-M 

20 


The run scripts are set to run through all 20 of the GCMs will the following line: 











































Handbook of Climate Change and Agroecosystems Downloaded from www.worldscientific.com 
by NASA GODDARD SPACE FLIGHT CENTER LIBRARY on 04/14/15. For personal use only. 


414 


Guide for Running AgMIP Climate Scenario Generation Tools With R 


Running AgMIP Cl ima te Scenario Generation scripts - 25 


run. gems <- 1:20 

The script is designed to run through the GCMs as they are ordered in the variable gemname (the 
same ordering as appears in the table above). The ordering of the GCMs is also provided in the 
commented lines directly above this input variable (as can be seen in the screenshot). Say, for 
example, you want to run the script for BNU-ESM2M, HadGEM2-CC and MIROC5 which are 
positioned as the 3 rd , 10 th and 15 111 in gemlist. You would then define as 

run. gems <- c (3, 10,15) 

If you only want to run the first five GCMs follow the same procedure and enter the line 

run. gems <- c(l,2,3,4,5) 

Or alternatively, 

run . gems <- 1:5 

run.rcps - This input variable, identified by the yellow arrow, is used to control which 
emissions scenarios (RCPs) will be used to create the output files. Below is a table denoting the 
values of run.rcps along with the corresponding emissions scenario (RCPs). 


RCP 

run.rcps value 

Historical 

1 

RCP 2.6 

2 

RCP 4.5 

3 

RCP 6.0 

4 

RCP 8.5 

5 


NOTE: While the script is designed to run RCP 2.6 (2) and RCP 6.0 (4), the data for these 
two emissions scenarios has not been included in the data provided by the l inks in the 
subsection Downloading R scripts and data from WebDrive . Thus, you should only 
define run.rcps as a 1, 3, or 5 (historical, RCP 4.S, or RCP 8.5) for the script to function 

cwr a etly i 

As can be seen in the screenshot, the run script “run_agnrip_simple_delta.R” is set to run RCP 
4.5 (3) and RCP 8.5 (5) with the following line: 
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run. reps <- c(3,5) 


If you wish to produce files from another emissions scenario simply redefine run.rcps. For 
example, if you only want to run the script from the emissions scenario RCP 8.5 (5), set the input 
variable to 5 using the following line: 

run.rcps <- 5 

If you want to run the script for the historical scenario (1) as well as RCP 4.5 (3) and RCP 8.5 (5) 
you can use the line 

run.rcps <- c(l,3,5) 


IMPORTANT: Some of the run scripts ( “run apmin simple. mandv.R” in narticnlart 
require a large amount of local memory. As such, we recommend that you run one 
emissions scenario at a time for “run aemip simple mandv.R " If the script breaks 
while you are running one of the AgMIP Climate scripts due to memory issues, consider 
reducing the number of emissions scenarios or number of future the time periods. After the 
run has successfully completed, reboot your computer (to reset your local memory) and run 
the other time periods and emissions scenarios you will require. 


run.decs - Another input variable, as identified by the green arrow, used to control the output 
files produced by the run scripts. This input variable is used to control which time periods will 
be used. Below is a table of the three preset future time periods. 


Future time period 

run.decs value 

First year 

Last year 

Near-term 

1 

2010 

2039 

Mid-Century 

2 

2040 

2069 

End-of-Century 

3 

2070 

2099 


Similar to and run.rcps, this variable is easily adjusted. Below is a line demonstrating how to 
produce output files for Near-term and End-of-Century future time scenarios 


run . decs 


<- c(l,3) 
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If you wish to produce files over just the Mid-Century time period, use the line provided below. 

run. decs <- 2 


datashort - This variable is used to define your sub-region for WorldClim files and is used 
exclusively by the script “ ran agmip farmclimate-R ” Your region is determined by the table 
below and you should enter your region into the input variable datashort. 


Datashort 

Region Name 

S. Lat 

N. Lat 

W. Lon 

E. Lon 

EAfrica 

Eastern Africa 

20“ S 

22“N 

20“ E 

53“E 

WAfrica 

Western Africa 

3“N 

25 “N 

18*W 

15°E 

SAfrica 

Southern Africa 

36“S 

8°S 

11“E 

51”E 

Pakistan 

Pakistan 

23 *N 

38°N 

60" E 

80“E 

IGB 

Indo-Gangetic Basin 

15“N 

38“N 

65"E 

97“E 

SIndiaSriLanka 

S. India and Sri Lanka 

5“N 

30“N 

70“ E 

90"E 


If, for example, your region is located at 5°N, 35° E, you would enter the input variable as 

datashort <- 'EAfrica' 


due to your location in the subregion of Eastern Africa. 

sitelat - Another input variable used exclusively with “ ran agmip farmcIimate.R ”. This input 
variable defines the series of site latitudes for which you are creating climate baseline files. This 
variable should be defined in terms of degrees North such that the values range from -90 °N 
(South Pole) to 90 “N (North Pole). 

For example, below is the line from the run script that will produce 3 output .AgMIP files for 
forms at latitudes -0.7 °N, -0.6 °N and -0.75 °N respectively. 

sitelat <- c(-00.70, -00.60, -00.75) 
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IMPORTANT: AgMIP Climate Scenario Generation scripts since version 2.0 of the Guide 
extract information from the header of the .AgMIP files, thus the coordinates of the seedfile 
should not be included in the datashort - This variable is used to define your sub-region for 
WorldClim files and is used exclusively by the script “run_agmip_fannclimate.R”. Your 
region is determined by the table below and you should enter your region into the input 
variable datashort. 


sitelon - Same as datashort - This variable is used to define your sub-region for WorldClim 
files and is used exclusively by the script “run_agmip_farmclimate.R”. Your region is 
determined by the table below and you should enter your region into the input variable 

datashort. 


Datashort 

Region Name 

S. Lat 

N. Lat 

W. Lon 

E. Lon 

EAfrica 

Eastern Africa 

20° S 

22 °N 

20°E 

53 °E 

WAfrica 

Western Africa 

3°N 

25°N 

18°W 

15°E 

SAfrica 

Southern Africa 

36°S 

8°S 

11 °E 

51“E 

Pakistan 

Pakistan 

23°N 

38°N 

60”E 

80°E 

IGB 

Indo-Gangctic Basin 

15°N 

38°N 

65 °E 

97°E 

SIndiaSriLanka 

S. India and Sri Lanka 

5°N 

30°N 

70°E 

90°E 


If, for example, your region is located at 5 ° N, 3 5 ° E, you would enter the input variable as 
datashort <- 'EAfrica' 
due to your location in the subregion of Eastern Africa. 


sitelat, but with longitude. This variable should be defined in terms of degrees East such that 
90° W should be entered either as -90 °E or as 270 °E. For example, below is the line from the 
run script. 


sitelon <- c(37.54, 37.58, 37.69) 


This example denotes that the run script will produce 3 output .AgMIP files for farms at 
longitudes 37.54 °E, 37.58 °E and 37.69 °E respectively. 


Modifying variables in the run scripts 
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NOTE: For the puiposes of this tutorial, we recommend that you perform the following 
modifications of the AgMIP Cl ima te scripts using the R GUI to gain some familiarity with 
R’s user interface. However, if you are familiar with R and feel more comfortable 
amending the scripts in your preferred text editor, feel free to do so. 


This subsection explains the steps that are required to prepare a run script for use on your 
computer. This subsection generalizes how to make the adjustments to the script The 
subsection. Suggested workflows for the run scripts , explains specifically how to run each of 
the AgMIP Climate run scripts. Like the previous subsection, Explanation of input variables 
to be changed in the run scripts . We reco mm end that you read through this subsection to get 
an idea of how to modify the run scripts and refer to it later if you have questions when working 
through the workflows presented in the subsection. Suggested workflows for the run scripts . 
The following section will explain: 

1 . How to open a script in R 

2. How to modify rootDir to tailor the run script for your computer 

3 . How to modify the other input variables 
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We will begin by opening the run script you are interested in running. After opening R, there are 
a couple ways you can open the run script You can either 

* Select “File > Open script. . from the menu bar along the top of the R window as 
indicated by the red arrow, or 

* Press the open button, in the toolbar (hidden in the below image by the drop-down 
window). 



If you followed the recommended folder structure outlined in the subsection. Creatine folders 
for R scripts and input data, the scripts will be located in the folder “~/R/r /”. 
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Once you have opened the run script in the R editor, the next step is to adjust rootDir to reflect 
the location of your “~/R/” folder. Scroll down section “Input variables” below the run script 
heading to the line that reads: 

rootDir <- ’*** your directory here ***\\R\\' 

Enter the path to your root directory, that is the location of your “~/R J" folder, discussed in the 
subsection. Creating folders for R scripts and input data in the run script using double 
forward slashes (\\) to separate the folders. For an example, refer to rootDir. 


NOTE: If you will be running these scripts for multiple station locations, we recommend 
that you save the run script at this point in time. The input variable rootDir should not 
change from run to run and thus you can save the run script to avoid having to adjust this 
variable for every run. The other input variables will be changed for different runs, so 
saving the script here is a logical saving point. 


After rootDir is assigned, you will not have to change any of the other folder locations as the 
folder structure is hardwired into the scripts. The AgM IP Climate scripts were intentionally 
designed this way to make this process as simple as possible. 

After saving the run script, you will have to define the input variables (described in the previous 
subsection. Explanation of input variables to be changed in the run scripts’) prior to every 
run. You can edit the script either in R editor or another text editor of your choosing. We will 
continue to use R editor in this Guide. 

Each AgM IP Climate Scenario Generation script has a section for the input variables that must 
be assigned to run the script successfully. The next section, Suggested workflows for the run 
scripts, outlines the input variables that will need to be defined. These input variables will need 
to be defined prior to every run. After properly defining these input variables, you would be 
ready to run, or source, the AgMIP climate script. 


IMPORTANT: At this point you have to save the run script to “~/R/r/”. We recommend 
that you save the script as a separate file with a unique name (e.g. 
“run_agmip_farmclimate_stnl” or “run_agmip_simple2full_USAM_delta3”) to reflect the 
particular station, dataset or analysis. This will retain the initial AgMIP run script as a 
template for future analyses. 
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Sourcing the ran scripts 


After defining the input variables of your run script and saving the run script with a unique name 
to u ~/RJrr, return to the R Console window. To run your amended run script, select “File > 
Serose R code... "from the menu bar as is indicated by the red arrow in the screenshot below. It 
is important that you are in the R Console, otherwise "Source R Code” will not be an option in 
the File dropdown menu. 


IMPORTANT: You must save the run script with your changes to the input variable 
definitions for the scripts to be sourced correctly. 



After you select “Source R code.. navigate to “~/R/r/” and select your antsided and saved run 
script Typically the run scripts will lake a few hours to run as they are batch creating multiple 
files based on a variety of scenarios. 
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The nest section. Explanation of outwit flics, describes the naming convention that is used with 
the 8 digit AgMIP file names. The following section. Suggested workflows for the AgMIP 
CHmate r un scripts, outline suggested procedures for using the different AgMIP run scripts. 

Explanation of output flies 

The provided AgMIP run scripts exeats a series of files that are named with an eight digit code. 
This naming convention was created to help convey details of the scenarios so that they may be 
understood by other AgMIP participants and the AgMIP IT tools. 

Below is a brief explanation of the key (pages 33-34) for this code. The file in the figure below is 
the first file that would be created from running “ran mrnitn lhnnle delta. H* 1 with the preset 
parameters using the given test AgMIP baseline file, ‘TJSMAOXXX. AgMIP”. 



The first 4 digits, as described by S digit baseline .AgMIP file name in the subsection Creating 
baseline AgMIP fflee with the enabled wacrn and haaeffla in the snWtinn Krai— Bon of 
Input variable* t» he dumped In the ran scripts, are used to describe the location of the 
climatological (meteorological) station where your historical dataset was collected. In the above 
example, the code “USAM” is used to describe the location where the first 2 letters are an 
abbreviation for the country (US = United States). The website, < 3ittp://www.web- 
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l.com/countrv-codes/ >. provides 2 letter codes for countries and we recommend that you use 
these abbreviations for the Erst two digits of your 8 digit code. 

The second 2 digits refer to the site location and are up to you to define. For the code 
“USAMCAXA", the “AM” is used to denote Ames, Iowa. Use your discretion when choosing 
these digits. These digits should help create a unique code for your file so try to avoid redundant 
digits. 

The fifth digit describes the time period and emissions scenario. Generally, the run scripts are 
designed to loop through multiple time periods and emissions scenarios when producing new 
files and this digit is automatically assigned to reflect the parameters imposed upon the dataset. 

In the example on the previous page, the fifth digit, C, denotes that the file uses RCP 4.5 over the 
future time period from 2010-2039. A baseline AgMIP file has a zero (0) as the fifth digit 
de fining that the file’s time period runs from 1980-2009 and does not have an imposed emissions 
scenario. 

The sixth digit typically denotes the imposed GCM scenario from CMIP5. For the purposes of 
many of the analyses that will be conducted in the future, the CMIP5 GCM outputs will be used 
unless otherwise mentioned. In the example, the sixth digit is an “A”, which tells us that this 
particular file had a GCM scenario imposed on the baseline data from ACCESS1-0. A baseline 
AgMIP file has an X as the sixth digit denoting that there was no GCM used in the creation of 
the file. Similar to the fifth digit, the sixth digit is typically assigned within the AgMIP run script 
loops. 

The seventh digit is used to impose additional downscaling or sensitivity tests on the baseline 
dataset. Generally this variable is unchanged when using the scripts described in this Guide 
unless you impose downscaling or scenario methodology independently. As a result, this digit is 
usually X, denoting that there has been no additional downscaling. 

The eighth digit tells us the type of scenario that has been imposed for this particular file. For 
the example code on the previous page, “USAMCAXA”, the eighth digit, “A”, tells us that there 
has been a mean change imposed a calculated from a GCM. This digit is typically assigned by 
the loop in the run script, but in when using the script “run_agmip_simple2full.R” you will need 
to identify this variable to specify files that you would like to be expanded to full files. The 
eighth digit is an “X” for the baseline AgMIP files, defining the file to be a dataset of daily 
observations. 

The next two pages (pages 33-34) list the key for the fifth, sixth, seventh and eighth digits of the 
AgMIP file naming convention. The first four digits are user defined and thus do not require a 
key. 
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Fifth Digit is time period and emissions scenario: 


0 = 1980-2009 baseline 

C = RCP4.5 

1 = A2-2005-2035 (Near-term) 

D = RCP6.0 

2 = Bl-2005-2035 (Near-term) 

E = RCP8.5 

3 = A2-2040-2069 (Mid-Century) 

F = RCP2.6 

4 = B 1-2040-2069 (Mid-Century) 

G = RCP4.5 

5 = A2-2070-2099 (End-of-Century) 

H = RCP6.0 

6 = B 1-2070-2099 (End-of-Century) 

I = RCP8.5 

S = sensitivity scenario 

J = RCP2.6 

A = observational time period (determined in 

K = RCP4.5 

file) 

L-RCP6.0 

B = RCP2.6 2010-2039 (Near-tenn) 

M = RCP8.5 


Sixth Digit is source of 
baseline data (if baseline 
scenario): 

X = no GCM used 

0 = imposed values (sensitivity 

tests) 

Q = Bias-corrected MERRA 
T — NASA POWER 
U = NARR 

V = Bias-corrected CFSR 
W~ MERRA 

Y = NCEP CFSR 

Z = NCEP/DoE Reanalysis-2 


Sixth Digit is GCM (if CMTP3 
scenario - these are not 
required for AgMIP full 
assessments): 

X = no GCM used 

0 = imposed values (sensitivity 
tests) 

A = bccr 
B = ccema cgcm3 
C = cnrai 
D - csiro 
E = gfdl 2.0 
F = gfdl2.1 
G = giss er 
H = inmem 3.0 

1 = ipsl cm4 

J = miroc3 2 medres 
K = miub echo g 
L = mpi echamS 
Nl — nun cgcm2 
N = near ccsxn3 
O = near pcml 
P = ukmo hadcm3 


2010-2039 (Near-term) 
2010-2039 (Near-term) 
2010-2039 (Near-tenn) 
2040-2069 (Mid-Century) 
2040-2069 (MidCentury) 
2040-2069 (MidCentury) 
2040-2069 (MidCentury) 
2070-2099 (End-ofCentury) 
2070-2099 (End-ofCentury) 
2070-2099 (End-ofCentury) 
2070-2099 (End-ofCentury) 


Sixth Digit is GCM (if CMIP5 
scenario): 

0 = imposed values (sensitivity 
tests) 

A = ACCESS 1-0 
B = bcc-csml-1 
C = BNU-ESM 
D — CamESM2 
E = CCSM4 
F = CESM1-BGC 
G — CSIRO-Mk3 -6-0 
H = GFDL-ESM2G 

1 = GFDL-ESM2M 
J = HadGEM2-CC 
K - HadGEM2-ES 
1. = imn em4 

M = IPSL-CM5A-LR 
N = IPSLCM5A-MR 
0-MIROC5 
P = MMOC-ESM 
Q = MPI-ESM-LR 
R = MPI-ESM-MR 
S = MRICGCM3 
T = NorESMl-M 
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Seventh Digit is downscaling/scenerio 
methodology: 

X = no additional downscaling 

0 = imposed values (sensitivity tests) 

1 = WRF 

2 = RegCM3 

3 = ecpc 

4 = hnn3 

5 = crcm 

6 = mm5i 

7 = RegCM4 
A = GIST 

B = M arkS IM 
C = WM2 

D = 1/8 degree BCSD 
E - 1/2 degree BCSD 
F - 2.5miniite WorldClim 
W = TRMM 3B42 
X = CMORPH 
Y-PERSIANN 
7. — GPCP 1DD 


Eighth Digit is Type of Scenario; 

X = Observations (no scenario) 

A = Mean Change from GCM 

B = Mean Ch ang e from RCM 

C - Mean Change from GCM modified by RCM 

D = Mean Temperature Changes Only 

E = Mean Precipitation Changes Only 

F = Mean and daily variability change for Tmax, 

T tnin andP 

G = P, Tmax and Train daily variability change 
only 

H = Tmax and T min daily variability and mean 
ch ang e only 

I = P daily variability and mean change only 
J = Tmax and Tmin daily variability change only 
K - P daily variability change only 
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Suggested workflows for the run scripts 

The following section details the specific workflows that should be used to run the AgMIP 
Climate Scenario Generation scripts. The previous section. Running AgMIP Climate Scenario 
Generation scripts, has generalized how to run these scripts while this section will present 
detailed workflows for the AgMIP Climate Scenario Generation scripts. 

If you have any questions about the input variables, refer to the subsection. Explanation of 
input variables to be changed in the run scripts or simply click on the input variable for which 
you have a question. 

In the overall workflow, this subsection picks up after you have successfully created your 
baseline .AgMIP file as described in the section Creating haseline .AglVTlP files with “AgMIP 
Excel File Template v2.0.xls” . This subsection then details how to: 

1. Create new baseline AgMIP files for sites in the region of your climate station using 
“run agmip farmclimate.R ”. 

2. Impose mean changes on the climate variables temperature maximum (Tmax), 
temperature minimum (Tmin) and precipitation (Rain) based on the GCMs with different 
emissions scenarios and across different time scales using 

“run agmip simple delta.R ". 

3. Impose mean and variability changes to the climate variables temperature maximum 
(Tmax), temperature minimum (Tmin) and precipitation (Rain) based on the GCMs with 
different emissions scenarios and across different time scales using 

“run agmip simple mandv.R ”. 

4. Convert file simple output files created by either “ run agmin simple delta.R ” or 
“run agmin simple mandv.R ” to full files (including wind speed (Wind) and relative 
humidity (Rhum) from the baseline AgMIP file and modified dew point (Dewp), vapor 
pressue (Vprs) derived from changes to Tmax) using “ run agmin simple2full.R ” 

To reiterate, the run scripts are designed to batch produce a number of . AgMIP files by calling 
the script of the same name but without the “run” prefix (for example, 

“ run agmip farmclimate.R’ ' calls “agmip farmclimate.R” and “run aemip simple2full.R ” 
calls “agmip_simple2full.R”). The tun scripts are designed to loop through multiple scenarios to 
reduce the number of times you will need to run these scripts. The following workflows explain 
how these run scripts should be used and have some information on how they can be tailored for 
your needs. 
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runagmipfarmclimate.R 

The run script “run agmip farmclimate.R ” creates baseline .AgMIP files for a series of 
locations (farms) in a given region. This script can be used to create climate files for a large 
number of farm locations that are close to a central weather station where the historical climate 
series has been quality controlled. If the farms are greater than ~100km away it may be better to 
pull the sites directly from Ag-MERRA dataset. 

The workflow is fairly simple if you have placed all of the files in the correct location as 
described in Creating folders for R scripts and input data . All you will have to do is switch a 
couple of the input variables in the run script prior to each run. 


To run this script, 

1 . Open R and then open the run script “ run agmip farmclimate.R ” to edit the input 
variables. You could also edit the run script in a text editor if you prefer. 

2. Change the rootDir to reflect the location of the folder “R” and datashort to reflect your 
WorldClim sub-region. 

3. Save the updated run script as your template “ run agmip farmclimateJt ” script. 

Adjust the input variables seedfile, shortregion, headerplus, datashort - This variable is used 
to define your sub-region for WorldClim files and is used exclusively by the script 
“run_agmip_farmclimate.R”. Your region is determined by the table below and you should 
enter your region into the input variable datashort. 


Datashort 

Region Name 

S. Lat 

N. Lat 

W. Lon 

E. Lon 

EAfrica 

Eastern Africa 

20° S 

22°N 

20°E 

53 °E 

WAfrica 

Western Africa 

3°N 

25°N 

18°W 

15°E 

SAfrica 

Southern Africa 

36°S 

8°S 

11°E 

51°E 

Pakistan 

Pakistan 

23 1 ’N 

38°N 

60°E 

80°E 

IGB 

Indo-Gangetic Basin 

15°N 

38°N 

65 "E 

97°E 

SIndiaSriLanka 

S. India and Sri Lanka 

5°N 

30°N 

70°E 

90°E 


If, for example, your region is located at 5 ’ N, 35 ° E, you would enter the input variable as 


datashort <- 'EAfrica' 


due to your location in the subregion of Eastern Africa. 
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4. sitelat, and sitelon to reflect your station information (of the base file). For headerplus, 
simply enter the location of the station (for example, “Ames, Iowa, USA” or “Nioro, 
Senegal). 

5. Save the amended script with a unique name (for example, 

“run_agmip_farmclimate_l.R” or “run_agmip_farmclimate_Stationl44.R”) to “~/R/r/”. 

6. Confirm the presence of your baseline file, “****OXXX. AgMIP” in 
“~/R/ data/Climate/Historical/”. 

7. Source (run) the script from the R Console. 

8. Try not to use your computer if possible to improve run time of script. 

9. Alter the script has completed, the baseline .AgMIP files for a series of locations should 
be located in “~\\R\\data\\Climate\\Historical\\”. Confirm the presence of these files. 

Although we did not provide you with a test file for this script, there is an example of an output 
from “ run agmin farm climate. R ” shown in Appendix 2 . The example file demonstrates the 
format of the .AgMIP files produced by this run script. This script should run quickly but will 
take more time depending on the number of sites for which you are producing baseline .AgMIP 
files. 


Note: When closing R, you will be asked to save your workspace. This saves the working 
environment (data, functions, etc.) and should not be necessary for the purposed of running 
the AgMIP Climate scripts. 


run_agmip_simpIe_delta.R 

This script creates delta scenarios from CMIP5 GCMs and BCSD in the AgMIP standard format 
You can either use your baseline AgMIP file or the newly created AgMIP files produced by 
“ rnn agmip farmcIimate.R ” as the input basefile for this run script. We are currently working 
on a run script that would loop through a series of basefiles to simplify larger batch runs. 

To run this script, 

1 . Open R and then open the run script “ run agmip simple delta.R” to edit the input 
variables. You could also edit the run script in a text editor if you prefer. 

2. Change the rootDir to reflect the location of the folder “~/R”. 

3. Save the updated run script as your template “ run agmip simple delta.R” script. 

Adjust the input variables basefile and 
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seedfile - Similar to basefile, seedfile is the 8 digit code of your baseline climate file. The input 
variable seedfile is used in the script “runagmipfarmclimate.R” as this will be the dataset 
from which other local farm sites’ climate files will be constructed. For example, the 8 digit 
code of the seedfile would be “KEMBOXXX” if the baseline AgMIP file was named 
“KEMBOXXX.AgMIP”. 

shortregion - Another input variable used exclusively with “mn_agiiiip_farmclimate.R” > 
shortregion is a 2 digit identifier used for a particular region for which you will be producing 
multiple files. Typically, it is sufficient to simply strip off the country code. For example, the 
country code “KE” can be stripped off the seedfile “KEMBOXXX” leaving the regional code 
“MB”. The resulting output files from “mn_agmip_farmcllmate.R” will be named 
“MBO 1 OXFX. AgMIP”, “MB020XFX.AgMIP”, “MB030XFX. AgMIP”. . . 

end.code - This input variable, used exclusively in “run_agmip_simple2full.R”, denotes the 
last two (seventh and eighth) digits of the files to be adjusted. For example, if you want to adjust 
all of the files of the structure “USAM**XA. AgMIP” with a baseline file of 
“USAMOXXX. AgMIP” to the full format, the 

would be “XA”. This would set the run script to cycle through the scenarios (defined by , 
run.rcps, and run.decs) of all the files in “~/R/data/Climate/Simplescenario” that begin with 
“USAM” and end with “XA”. 


basedecs to reflect the relevant information of your baseline AgMIP file. You can also 
adjust the variables , run.rcps, and run.decs if you only require a subset of these input 
variables. 


IMPORTANT: If you are having issues completing the script due to local memory issues, 
consider adjusting , run.rcps, and run.decs. Running subsets will produce fewer output 
files and use less of your local memory. After completing such a subset it is important to 
reboot your computer prior to running the next subset. 


4. Save the amended script with a unique name (e.g. “run_agmip_simple_delta_Stnl2.R” or 
“run_agmip_simple_delta_Farm67.R”) to “~/R/r f\ 

5 . Confirm the presence of your basefile in “~\\R\\data\\Climate\\Historical\\” 

6. Source (run) the script from the R Console 

7. Try not to use your computer if possible to improve run time of script. 
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8. After the script has completed, the basic delta adjusted files should be located in 
~\\R\\data\\Climate\\Simplescenario\\. Confirm the presence of these files. 

The .AgMIP files produced by “ run amnio simple delta.R” only contain the climate variables 
Srad, Tmax. Tmin, and Rain. If you require any of the other climate variables (Wind, Dewp, 
Vprs, and/or Rhum) you can use “run agmip simple2full.R ” to create full files. The 
subsection Creating “full” files from “run apmin simnle delta Jt” outputs explains how to 
do this. 

If you chose to source “run agmip simple delta.R ” with the test file, “USAMOXXX.AgMIP” 
to produce the simple future scenario files in .AgMIP format, you can check the first created file, 
“US AMCAXA.AgMEP”, against the data shown in Appendix 3 . However, the data shown in 
this Appendix is the output from “ run agmin simnle delta Jt ” and 
“ run agmin simnle2fullJt” so remember that the example displays data for wind speed 
(WIND), dew point (DEWP), vapor pressure (VPRS), and relative humidity (RHUM) as well. 


run_agmip_simple_mandv.R 

The script “agmip simple mandv.R” applies mean and variability changes to temperature 
maximum (Tmax), temperature minimum (Tmin) and precipitation (Rain) of the input baseline 
.AgMIP file. The other climate variables (Srad, Wind, Dewp, Vprs, Rhum) are not adjusted by 
this script. 

Running this script is similar to “run agmip simple deltaJR ” in that you should only have to 
change the root directory (rootDir), adjust the input variables, and the run script will impose a 
series of changes to the baseline file to create a batch of output files. Below is a suggested 
workflow for running “ ran agmin simnle mandv.R ”: 

To run this script, 

1 . Open R and then open the run script “ ran agmip simple mandv.R ” to edit the input 
variables. You could also edit the run script in a text editor if you prefer. 

2. Change the rootDir to reflect the location of the folder “~/R”. 

3. Save the updated run script as your template “ run agmin simnle delta.R ” script. 

Adjust the input variables basefile and headerplus to reflect the relevant information of 
your baseline AgMIP file. Change headerplus from “Ames, Iowa, USA” to the text with 
the station or site location. We recommend that you leave the text “with mean and daily 
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variability changes for Tmax, Tmin and P”. You can also adjust the variables , run.rcps, 
and run.decs to run a subset of these input variables. 


IMPORTANT: This run script, “ run agmip simple mandv.R ’' will quickly use up your 
local memory. We recommend that you only run one emissions scenario (as defined by 
run.rcps) and one time period (as defined by run.decs) per run to ensure the script finishes 
without encountering any memory issues. It is important to remember to restart your 
computer after completing a run if you intend on running another analysis. 


4. Save the amended script with a unique name (e.g. “run_agmip_simple_mandv_Stnl2.R” 
or “run_agmip_simple_mandv_Farm67.R") to “~/R/r/”. 

5 . Confirm the presence of your basefile in “~\\R\\data\\Climate\\Historical\\” 

6. Source (run) the script from the R Console 

7. Try not to use your computer if possible to improve run time of script. 

8. After the script has completed, the basic delta adjusted files should be located in 
~\\R\\data\\Climate\\Simplescenario\\. Confirm the presence of these files. 

The . AgMIP files produced by “run agmin simple mandv.R ” only adjust the climate variables 
Tmax. Tmin, and Rain. The other climate variables (Srad, Wind, Dewp, Vprs, and Rhum) are 
pasted into this file from the AgMIP file you specified as the basefile. If you require that Dewp 
and Vprs are updated to reflect changes in Tmax, you can use “ run aymip simple2fnIl-R ” to 
create full files. The subsection Creating “full” files from “run agmip simple mandv.R r ’ 
outputs explains how to do this. 


IMPORTANT: The script “ agmin simple mandv.R” occasionally has difficulties 
applying the mean and variability changes without significantly changing the distribution and 
will notify you that there is a poor fit. The script is designed to then loop through and try to 
apply the changes two more times. If the run script is unable of applying the changes after 
three tries, the script will fill in that particular month with -99 values. 

If you encounter this issue, you can try to run the script again (after it has completed creating 
the other output files and rebooting your computer) specifically for this file by using , 
run.rcps, and run.decs. If the problem persists, please contact us so we can help you resolve 
this issue and produce the . AgMIP files you require for your research. 
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run_agmip_simple2full.R 

The script “agmip_simple2full.R” converts basic future scenarios (Srad, Tmax, T min, Rain) into 
full scenarios with relative humidity-controlled dew point and vapor pressure based upon daily 
Tmax. The script fills Wind and Rhum with the data from the input baseline AgMIP file and 
computes new values for Dewp and Vprs based on Tmax and Rhum. 

The run script can be used for output AgMIP files from “ run agmin simple delta.R ” and from 
“ run agmin simple mandv.R *’. Both have slightly different workflows, thus there are two 
subsections. Creating “full” files from “run agmin simple delta.R” outputs and Creating 
“fall” files from “run ap-mip simple mandv.R” outputs , describing how to use this run 
script. 


Creating “fall” files from “run agmip simple delta.R” outputs 

After running “run agmin simple delta.R ” you will have a series of basic delta adjusted 
.AgMIP files located in “~/R/data/Climate/Simplescenario/”. These new AgMIP files contain 
adjusted data for temperature maximum (Tmax), temperature minimum (Tmin) and precipitation 
(Rain) as well as the baseline data for solar radiation (Srad). 

If you require the full delta adjusted .AgMIP files with wind speed (Wind), dew point (Dewp), 
vapor pressure (Vprs), and the relative humidity (Rhum), you will have to run the script 
“ run agmin simnle2fall.R” This run script will paste the baseline values for Wind and Rhum 
and will calculate new values for Dewp and Vprs using relative humidity controlled vapor 
pressure based on daily Tmax.. 


IMPORTANT: The script, “run agmip simple2full.R ". uses the outputs from 
“run agmin simple delta.R ” to create the full delta adjusted .AgMIP files. Thus 
“run agmin simple delta.R ” must be run prior to “ run agmin simnle2fall.R ” for the 
script to function correctly. 


To run this script, 

1 . Open R and then open the run script “ run agmin simple2fulljt” to edit the input 


variables. You could also edit the run script in a text editor if you prefer. 
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2. Change the rootDir to reflect the location of the folder “~/R”. 

3. Save the updated run script as your template “ run aemip simple2full.R” script. 

Adjust the input variables baseflle and 

to reflect the relevant information of your baseline AgMIP file. For AgMIP files produced by 
“ run agmip simple delta.R” . the variable 

should be defined as “XA” unless you have made other downscaling adjustments to the 
files. You should also adjust the variables , run.rcps, and run.decs to run a subset of 
these input variables based on the files you produced by “ run aemip simple delta.R ” 
and the files for which you require a full climate dataset. 

4. Save the amended script with a unique name (e.g. “run_agmip_simple2fidl_l .R” or 
“run_agmip_simple2full_deltas.R”) to “~/R/r/”. 

5. Confirm the presence of your baseflle in “~/R/data/CIimatc/Historical/” as well as the 
simple files to be converted in “~/R/data/Climate/Simplescenario/”. 

6. Source (run) the script from the R Console 

7. Try not to use your computer if possible to improve run time of script. 

8. After the script has completed, the basic delta adjusted files should be located in 
~\\R\\data\\Climate\\Fullscenario\\. Confirm the presence of these files. 

If you chose to run “run agmlp simple delta.R ” and “ run aymip simple2f»ll.R ” with the 
test file, “USAMOXXX. AgMIP” to produce the full future scenario files in .AgMIP format, you 
can check the first created file “USAMCAXA. AgMIP”, against the data shown in Appendix 3 . 

Creating “full” files from “run agmip simple mandv.R” outputs 

If you require that the mean and variability change is also applied dew point (Dewp) and vapor 
pressure (Vprs), the run script “ run aemin simnle2fnll.R ” can be used to expand these files. 
The script “ run agmip simnle2full.R ” was updated for version 2.0 of the Guide to use the files 
produced by “ run aemip simple mandv.R ” (found in “~/R/data/Clixnate/Simplescenario”). 
The run script calculates new daily values for Dewp, and Vprs based on the baseline relative 
humidity (Rhum) and the newly calculated temperature maximum (Tmax). The run script 
produces new “full” files located in “~/R/data/Climate/Fullscenario”. 

To run this script and convert the simple files to full, 

1 . Open R and then open the run script “ run aemip simpIe2full.R” to edit the input 
variables. You could also edit the run script in a text editor if you prefer. 

2. Change the rootDir to reflect the location of the folder “~/R”. 
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3. Save the updated run script as your template “ run acnmin simnlc2full.R ” script. 

Adjust the input variables basefile, and 

to reflect the relevant information of your baseline AgMIP file. For AgMIP files produced by 
“ run agmin simnle mandv.R ”. the variable 

should be defined as “XF” where “X” denotes “no additional downscaling” and “F’ denotes 
“Mean and daily variability change for Tmax, Tmin, and P”. If you have made additional 
downscaling adjustments to your AgMIP files to be converted, adjust 

accordingly. You can also adjust , run.rcps, and run.decs to run a subset of these input 
variables based on the files you produced by “ run agmip simple mandv.R ” and the 
files for which you require a full climate dataset. 

4. Also, you should change the headerplus information to read “<your site location, 
country> with mean and daily variability changes for Tmax, Tmin and P” 

5. Save the amended script with a unique name (e.g. “run agmip simplc2full_2.R” or 
“ran_agmip_simple2full_mandv_l.R”) to “~/R/r/”. 

6. Confirm the presence of your basefile in “~/R/data/Climate/Historical/” as well as the 
simple files to be converted in “~/R/data/Climate/Simplescenario/”. 

7. Source (run) the script from the R Console 

8 . Try not to use your computer if possible to improve ran time of script. 

9. After the script has completed, the basic delta adjusted files should be located in 
~\\R\\data\\Climate\\Fullscenario\\. Confirm the presence of these files. 
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Next Steps... 

This concludes Version 2.0 of the Guide for R unnin g AgMIP-Climate Scenario Generation 
Tools with R in Windows. We will be expanding this Guide in the future to include other climate 
scenario tools. Your feedback is helpful and we appreciate any suggestions or recommendations 
you may have to improve this Guide. 

If you have further questions, contact please feel free to contact Nicholas Hudson at 
nih21()6(S).columbia.edn . 

T hanks again for your collaboration and good luck with R! 
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Appendix 1: Creating .AgMIP files in Excel 

This appendix describes how to create .AgMIP files in Excel manually without use of “AgMIP 
Excel File Template_v2.0.xls”. 

Here's the general approach: 

1) Arrange data series in excel spreadsheet with propter column order, proper units, etc. Missing 
values should be filled with -99 (you can get RHUM and wind speed from the MERRA-based 
dataset I provided if you do not have these, but Tdew and VPRS will need additional conversion 
and can be left as missing for the purposes of the R routines described in this Guide. 

2) Highlight all of the four columns that indicate date: @DATE, YYYY, MM, and DD (e.g., 
their first entries are "1980001 1980 1 1"), right click, and select "format cells". Select the 
number category and then change the number of Decimal places to 0. 

3) Highlight the columns corresponding to SRAD, TMAX, TMIN, RAIN, WIND, DEWP, and 
VPRS, right click, and select "format cells". Select the number category and then change the 
number of Decimal places to 1 . 

4) Highlight the column corresponding to RHUM, right click, and select "format cells". Select 
the number category and then change the number of Decimal places to 0. 

5) Right-click on the column letter at the top of the excel spreadsheet for the @DATE column 
(e.g., the "A" for Column A), select "Column Width" and set the column width to 7. 

6) Right-click on the column letter at the top of the excel spreadsheet for the YYYY column 
(e.g., the "B" for Column B), select "Column width", and set the column width to 6. 

7) Right-click on the column letter at the top of the excel spreadsheet for the MM and DD 
columns ("C" and "D"), select "Column width", and set the column width to 4. 

8) Right-click on the column letter at the top of the excel spreadsheet for the remaining columns 
(SRAD to RHUM; "E" to "L"), select "Column width", and set the column width to 6. 

9) Select "Save as" from the file menu and save as type "Formatted Text (Space 
Delimited)". Type in a full name with the extension (e.g., "USNYOXXX.AgMIP"). 
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10) You may need to click yes on a dialog box asking if you want to save only this sheet, and 
after you click that you can cancel any other dialog boxes or separate save as screens that may 
appear, 

11) Open the .AgMIP file in notepad and then copy 5-line header from a different existing 
.AgMIP file and paste it on top of the AgMIP file you just created. Modify header as necessary 
(e.g., header, 4-letter code, latitude, longitude, elevation), TAV is the average of all Tmax and 
Ttnin in the dataset, and AMP is half of the difference between the warmest average monthly 
temperature and the coolest average monthly temperature averaged over the entire dataset. 
REFHT is the height of the thermometer and WNDHT is the height of the anemometer (put -99 
if unknown). 

12) Save the file to record these modifications. 
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Appendix 2: Example file for “agmip farmclimate.R” 


Below is an example output data file created by “ run agmip farmclimate.R ” and located in 
“~\\R\\data\\Climate\\Historical\\”. This data file is provided so that you can confirm that your 
script ran correctly. The first 32 days of the data have been provided below as well as the last 5 
days of the data set (December 27-31, 2010). 


*’ WEATHER DATA : Embu, Kenya, cast to she 01 using WorldClim-derived climatological differences 


@ INSI 

LAT 

LONG 

ELEV 

TAV 

AMP 

REFHT 

WNDHT 





MB 01 

-0.700 

37.540 

1097 

20.4 

2.0 

-99.0 

-99.0 





@DATE 

YYYY 

MM 

DD 

SRAD 

TMAX 

TMIN 

RAIN 

WIND 

DEWP 

VPRS 

RHUM 

1980001 

1980 

1 

1 

23.0 

24.4 

12.6 

0.0 

238.4 

18.0 

20.9 

67 

1980002 

1980 

1 

2 

23.8 

27.1 

12.1 

0.0 

202.7 

19.6 

23.1 

63 

1980003 

1980 

1 

3 

21.5 

25.3 

10.1 

0.0 

234.7 

16.0 

18.3 

56 

1980004 

1980 

1 

4 

23.9 

26.2 

12.3 

0.0 

220.0 

15.9 

18.2 

52 

1980005 

1980 

1 

5 

24.3 

29.7 

13.7 

0.0 

197.3 

21.8 

26.5 

62 

1980006 

1980 

1 

6 

23.1 

25.9 

12.1 

0.0 

163.4 

17.7 

20.4 

60 

1980007 

1980 

1 

7 

23.2 

27.6 

10.3 

0.0 

173.3 

18.7 

21.8 

58 

1980008 

1980 

1 

8 

24.4 

27.6 

12.9 

0.0 

213.5 

16.9 

19.4 

51 

1980009 

1980 

1 

9 

21.8 

25.6 

12.9 

0.0 

181.7 

16.2 

18.6 

56 

1980010 

1980 

1 

10 

23.8 

26.7 

12.5 

0.0 

104.6 

14.7 

16.9 

47 

1980011 

1980 

1 

11 

23.5 

28.2 

12.2 

0.0 

186.3 

16.5 

18.9 

48 

1980012 

1980 

1 

12 

19.9 

25.9 

14.1 

0.0 

199.3 

15.3 

17.5 

51 

1980013 

1980 

1 

13 

22.7 

25.9 

13.1 

0.0 

196.4 

15.6 

17.9 

52 

1980014 

1980 

1 

14 

22.4 

26.2 

11.6 

0.0 

212.4 

16.5 

18.9 

54 

1980015 

1980 

1 

15 

22.4 

27.5 

12.7 

0.0 

254.4 

17.4 

20.1 

54 

1980016 

1980 

1 

16 

20.2 

25.9 

15.1 

0.0 

179.7 

17.1 

19.7 

58 

1980017 

1980 

1 

17 

23.6 

25.9 

12.6 

0.0 

255.5 

16.8 

19.3 

57 

1980018 

1980 

1 

18 

21.3 

27.1 

13.3 

0.0 

183.0 

16.8 

19.2 

52 

1980019 

1980 

1 

19 

23.2 

27 2 

13.6 

0.0 

187.3 

15.6 

17.8 

48 

1980020 

1980 

1 

20 

24.1 

27.8 

11.6 

0.0 

237.0 

17.1 

19.7 

51 

1980021 

1980 

1 

21 

23.3 

28.1 

12.7 

0.0 

281.5 

15.3 

17.6 

45 

1980022 

1980 

1 

22 

23.6 

27.1 

13.1 

0.0 

247.3 

16.1 

18.5 

50 

1980023 

1980 

1 

23 

22.7 

27.6 

14.6 

0.0 

188.2 

13.8 

15.9 

42 

1980024 

1980 

1 

24 

21.0 

29.5 

14.5 

0.0 

146.4 

16.6 

19.1 

45 

1980025 

1980 

1 

25 

22.2 

28.3 

15.4 

0.7 

107.6 

17.2 

19.9 

50 

1980026 

1980 

1 

26 

22.1 

25.2 

18.1 

0.3 

118.5 

16.7 

19.2 

59 

1980027 

1980 

1 

27 

19.5 

25.5 

15.0 

0.0 

155.6 

20.9 

25.1 

76 

1980028 

1980 

1 

28 

3.5 

25.9 

16.1 

0.0 

148.9 

19.7 

23.3 

68 

1980029 

1980 

1 

29 

18.4 

26.5 

14.6 

0.0 

136.2 

14.9 

17.1 

48 

1980030 

1980 

1 

30 

3.1 

25.8 

18.1 

0.0 

102.9 

15.2 

17.4 

51 

1980031 

1980 

1 

31 

2.9 

26.8 

16.7 

5.4 

127.4 

16.8 

19.3 

54 

1980032 

1980 

2 

1 

20.0 

26.0 

16.1 

0.0 

1504 

17.5 

20.2 

59 

2010361 

2010 

12 

27 

18.7 

26.0 

17.3 

5.8 

115 

16.2 

18.5 

54 

2010362 

2010 

12 

28 

24.0 

24.2 

16.7 

2.5 

71.1 

15.8 

18.1 

59 

2010363 

2010 

12 

29 

20.1 

26.0 

15.0 

0.0 

75.4 

14.4 

16.5 

48 

2010364 

2010 

12 

30 

19.8 

25.6 

13.9 

0.0 

78.8 

15.5 

17.7 

53 

2010365 

2010 

12 

31 

23.7 

26.2 

14.3 

0.0 

70.7 

13.2 

15.3 

44 
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Appendix 3: Example file for “agmipsimpledelta.R” and 
“agmip_simple2fuIl.R” 


Below is the data from “USAMCAXA.AgMIP”, the first file that is created by 
“run agroin simnle delta.R ” and “run anmin simnle2fiill-R ” located in 
“~\\R\\data\\Climate\\Fullscenario\\”. This data file is provided so that you can confirm that your 
script ran correctly. The first 32 days of the data have been provided below as well as the last 5 
days of the data set (December 27-31, 2010). 

’WEATHER DATA : USAMCAXA - baseline dates maintained for Leap year consistency 


@ INSI 

LAT 

LONG 

ELEV 

TAV 

AMP 

REFHT 

WNDHT 





USAM 

42.017 

-93.750 

329 

11.2 

14.6 

-99.0 

-99.0 





@DATE 

YYYY 

MM 

DD 

SRAD 

TMAX 

IMIN 

RAIN 

WIND 

DEWP 

VPRS 

RHUM 

1980001 

1980 

1 

1 

13 

13 

-1.5 

0.0 

3.1 

-03 

6.0 

89 

1980002 

1980 

1 

2 

4.7 

-0.3 

-2.6 

0.0 

4.9 

-7.6 

3.5 

58 

1980003 

1980 

1 

3 

1.9 

-0.3 

-4.8 

0.0 

4.3 

-9.0 

3.1 

52 

1980004 

1980 

1 

4 

3.8 

0 2 

-2.6 

0.0 

4.1 

-52 

43 

67 

1980005 

1980 

1 

5 

1.0 

02 

-3.2 

1.5 

3.4 

-25 

5.1 

82 

1980006 

1980 

1 

6 

8.5 

1.9 

-7.0 

2.1 

9.1 

-0.8 

5.7 

82 

1980007 

1980 

1 

7 

6.7 

-6.4 

-15.9 

0.0 

8.2 

-134 

23 

58 

1980008 

1980 

1 

8 

6.7 

-8.7 

-16.5 

0.0 

2.9 

-17.1 

1.6 

51 

1980009 

1980 

1 

9 

23 

-11.4 

-20.4 

1.5 

3.7 

-163 

1.8 

68 

1980010 

1980 

1 

10 

8.0 

63 

-13.7 

0.0 

8.0 

5.3 

8.9 

93 

1980011 

1980 

1 

11 

4.0 

11.3 

-9.8 

0.0 

11.9 

73 

103 

76 

1980012 

1980 

1 

12 

8.6 

0.8 

-13.7 

0.0 

6.1 

-10.4 

2.8 

43 

1980013 

1980 

1 

13 

8.6 

12.5 

13 

0.0 

6.7 

2.4 

73 

50 

1980014 

1980 

1 

14 

2.1 

6.9 

-3.2 

0.0 

4.1 

-1.0 

5.7 

57 

1980015 

1980 

1 

15 

1.0 

8.6 

4.6 

0.0 

4.1 

7.2 

10.2 

91 

1980016 

1980 

1 

16 

1.9 

7.5 

2.4 

24.4 

5.6 

7.4 

103 

99 

1980017 

1980 

1 

17 

3.7 

2.5 

-2.6 

1.2 

5.6 

-0.6 

5.9 

80 

1980018 

1980 

1 

18 

1.4 

13 

-3.2 

0.0 

2.7 

1.2 

6.6 

99 

1980019 

1980 

1 

19 

1.3 

13 

-4.3 

2.1 

43 

-0.0 

6.1 

91 

1980020 

1980 

1 

20 

3.1 

13 

-8.2 

0.0 

2.1 

-6.1 

3.9 

58 

1980021 

1980 

1 

21 

10.4 

13 

-7.0 

0.0 

4.2 

-0.0 

6.1 

91 

1980022 

1980 

1 

22 

8.8 

1.9 

-5.4 

0.0 

8.2 

-13 

5.6 

80 

1980023 

1980 

1 

23 

2.9 

-4.2 

-10.9 

0.0 

4.0 

-13.4 

23 

49 

1980024 

1980 

1 

24 

5.0 

6.9 

-7.0 

0.0 

5.7 

3.9 

8.1 

81 

1980025 

1980 

1 

25 

6.9 

5.8 

-4.8 

0.0 

5.9 

3.5 

7.8 

85 

1980026 

1980 

1 

26 

7.5 

-4.2 

-143 

0.0 

6.1 

-11.5 

2.6 

57 

1980027 

1980 

1 

27 

5.9 

-9.8 

-13.2 

0.0 

4.5 

-163 

1.8 

60 

1980028 

1980 

1 

28 

9.0 

-8.7 

-12.6 

0.0 

4.6 

-14.7 

2.0 

62 

1980029 

1980 

1 

29 

4.8 

-9.2 

-16.5 

0.0 

3.8 

-14.6 

2.0 

65 

1980030 

1980 

1 

30 

12.7 

-9.8 

-12.6 

1.5 

4.7 

-12.6 

23 

80 

1980031 

1980 

1 

31 

4.7 

-10.9 

-19.8 

0.0 

3.8 

-133 

23 

83 

1980032 

1980 

2 

1 

6.7 

-8.5 

-22.8 

0.0 

3.0 

-11.4 

2.6 

80 

2010361 

2010 

12 

27 

3.5 

1.1 

-0.3 

0.0 

2.0 

-0.8 

5.8 

87 

2010362 

2010 

12 

28 

4.8 

1.7 

-0.1 

9.0 

1.4 

-0.7 

5.8 

84 

2010363 

2010 

12 

29 

3.3 

2.8 

13 

2.8 

4.2 

1,8 

7.0 

93 

2010364 

2010 

12 

30 

2.6 

4.0 

1.7 

0.0 

5.1 

2.0 

7.1 

87 

2010365 

2010 

12 

31 

3.4 

13 

-1.5 

0.9 

5.3 

-25 

5.1 

76 
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